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The low and high-temperature phases of V4O7 have been studied by ab initio calculations. At 
high temperature, all V atoms are electronically equivalent and the material is metallic. Charge 
and orbital ordering, associated with the distortions in the V pseudo-rutile chains, occur below the 
metal-insulator transition. Orbital ordering in the low-temperature phase, different in V'^"'" and 
V^^ chains, allows to explain the distortion pattern in the insulating phase of V4O7. The in-chain 
magnetic couplings in the low-temperature phase turn out to be antiferromagnetic, but very different 
in the various V*"^ and V^"^ bonds. The V*"*" dimers formed below the transition temperature form 
spin singlets, but V'^^ ions, despite dimerization, apparently participate in magnetic ordering. 

PACS numbers: 71.30.-|-h, 71.20.-b, 75.25.Dk 



I. BACKGROUND 

The transition from a localized to an itiner- 
ant electron picture has been thoroughly stud- 
ied in the last decadeSfl In particular, 3d elec- 
tron systems can be described by both pictures, 
depending on the situation, structure, U/W ra- 
tio (on-site Coulomb repulsion vs. bandwidth), 
pressure, strain, etcj^ Metal-insulator transitions 
(MIT) have been observed in these compounds in- 
duced by pressure)^ chemical pressure;^ change in 
structure,— temperature^ spatial confinement,— or a 
simple change in electron count .--JJ 

Charge ordering (CO) plays an important role in 
the insulating side of these MIT in various tran- 
sition metal oxides, happening both in bulk com- 
pounds such as perovskite manganites^ii and nick- 
elatesi^, and in nanostructuresF^^ When materials 
approach a MIT, the lattice can become softer due 
to bond- length fluctuations,^'* and this could lead to 
structural distortions or phase transitions of peculiar 
character li^ii^ 

V4O7 is a member of the homologous series 
V„02„-i= V2O3+ (n-2)V02 (3< n <9), known as 
Magneli phase compounds The structure of the 
members of this series is based on the rutile unit 
cell. VO2, like several other transition metal diox- 
ides, crystallyzes (above room temperature) in the 
rutile structurCfi^ characterized by an octahedral 
coordination of the V cations by O atoms and a 
strong metal-metal bonding along the c-direction. 
The crystal structure of V4O7 (Ref. 19) consists of 
rutile-like blocks that are infinite in two-dimensions 
and four VOe octahedra wide in the third (see Fig. 



[T]). The shear planes between the rutile blocks have 
the typical local structure for V2O3 (Ref. [13) • The 
resulting structure is triclinic, belonging to a PI 
space group, with four and seven independent V 
and O sites, respectively. The V sites are split into 
two groups: V3 and V4, which are at the shear 
plane, and VI and V2, at the center of the rutile 
blocks. Since V3 and V4 octahedra are connected 
by face sharing, their environment is sesquioxide-like 
whereas that of VI and V2 is rutile-like. The cations 
form two independent chains V3-V1-V1-V3 and V4- 
V2-V2-V4 running parallel to the pseudorutile c-axis 
(see Fig. [1]). 

Various members of the series are characterized 
by a MIT.- The physics of the transition in VO2 
has been strongly debated. This is a correlated 
compound, but the occurrence of a dimerization 
along the c-direction, leading to a formation of spin- 
singlets due to the Peierls distortion, shows that 
both correlations and Peierls physics play a role in 
causing the transitioni^i"— V4O7 shows a MIT at 
250 K ijc)^ The MIT in V4O7 is a weakly first 
order onci^^^ In this case, the lattice symmetry 
is unchanged during the phase transition unlike the 
cases of V2O3 and V02. The compound also exhibits 
a paramagnetic to antiferromagnetic (AFM) transi- 
tion at Tat- 40 K.^^ The fact that the MIT occurs 
at a different temperature than the magnetic one al- 
lows to think that the magnetic order is not a driv- 
ing force of the phase transition. A simple electron 
count in V4O7 gives an average valence V"^'^+ for the 
metal cations. Hence, 2V'^+(d^) and 2V''+(d*) per 
formula unit could be expected if CO would occur. 
Different previous works showed that the metallic 



state was characterized by an almost complete disor- 
der of the and V''"'" cations whereas the insulat- 
ing state showed a high degree of ordering!^"— the 
V3-V1-V1-V3 chain is formed by ions and the 
V4-V2-V2-V4 by ones. To summarize, CO in 
the cation sites occurs below the transition, and this 
CO state eventually orders magnetically at lower 
temperatures. Other works point out that the mag- 
netic moments are localized not only in the insulat- 
ing phase but to some degree also in the metallic 
one.^^ Most studies show that, in addition to elec- 
tron localization, half of the V''"'" and all V'^"'" cations 
are paired, forming short V-V bonds along their cor- 
responding chains in the insulating state J^^^— An- 
other feature of the structural changes occuring at 
the MIT is the fcrro- or antiferroelectric-like shift of 
unpaired V^"^ ions at the shear planes, from the cen- 
ter of Og octahedra towards one of the surrounding 
oxygens. Such distortion ("twisting") also occurs in 

V02^'^"-2^ 

Several questions arise when looking at these re- 
sults. The first is, what is the nature of V pairing 
in the insulating phase? Similar pairing in VO2 is 
largely connected with the orbital repopulation)^ 
d-electrons, which in the metallic phase above Tc 
were more or less equally distributed over all the t2g- 
orbitals, below Tc occupy predominantly one orbital 
{dxy in the local coordinate system), with the lobes 
in the c-direction, giving strong overlap with neigh- 
boring V's along the chains in the c-axis. The re- 
sulting one-dimensional electron system is unstable 
with respect to Peierls, or spin-Peierls dimerization, 
leading to the formation of singlet V''"'"-V^"'" dimers. 
There exist many other compounds with transition 
metal ions with configuration d"'^ (containing V^^, 
Ti'^+ , or Nb''+ ) which also form singlet dimers at low 
temperatures. Thus, we could think that a similar 
phenomenon could exist also in V4O7, at least for 
the chains V3-V1-V1-V3. But, in fact, in these 
segments only the two middle V's (Vl-Vl) form sin- 
glet dimers, with V3 ions remaining unpaired and 
apparently magnetic. The question is why is this so. 

On the other hand, all the V'^"'" in V4O7 form 
dimers, which is not very typical for V^"*" cations 
(d^, S=l): only very few compounds with V'^+ form 
dimers. This happens for example in ZnV2 04,'^* but 
in this case the short V-V dimers are not in a singlet 
state, but are "ferromagnetic" (FM), i.e. have par- 
allel spins on two V's. What is the character of V^"*" 
dimers in V4O7, and why these dimers are formed, 
is an open question. 

Another problem is the character and the mecha- 
nism of magnetic ordering in V4O7 below Tjv = 40 
K. Despite long history of study of V4O7, its mag- 
netic structure is actually not known. It is also not 
clear which V ions participate in the magnetic or- 



dering. Apparently, the paired V*+ ions (Vl-Vl) 
are in a singlet state and "drop out of the game" . 
If V'^+ pairs would also be singlet, then only 1/4 of 
all the V's, the unpaired V^"*" in V3 position, would 
be magnetic. But then it would be difficult to un- 
derstand why such a "dilute" magnetic system, with 
magnetic ions with small spin S= 1/2 and located 
far from each other, would give a reasonably high 
value of the Neel temperature of 40 K. Alterna- 
tively, it is possible that the spins of the V'^'*' ions 
are not completely quenched, and that these ions 
also contribute to magnetic ordering. 

In order to shed light on the issues of charge and 
magnetic order and on the nature of the observed 
structural distortions, to analyze what degree of 
charge disproportionation exists in both phases and 
also what is the nature and strength of each of the 
magnetic couplings, we have studied the electronic 
structure of both metallic and insulating phases of 
V4O7 by means of ab initio calculations)^ Besides 
CO, we have obtained the specific orbital ordering 
(00) in the insulating phase, which allows to ex- 
plain the distortion pattern observed in it. We have 
also analyzed the plausible magnetic configurations 
and have discussed the magnetic properties of both 
low and high-temperature phases. 



II. COMPUTATIONAL DETAILS 

Our electronic structure calculations have been 
performed with the WIEN2K coder^^ based on 
density functional theory^i*^ (DFT) utilizing the 
augmented plane wave plus local orbitals method 
(APW-hlo)42. We have used the LDA-I-U'*^ approach 
including self-interaction corrections in the so-called 
"fully localized limit" with the on-site Coulomb re- 
pulsion U= 4.1 eV and the on-site Hund's rule cou- 
pling J= 0.7 eV. This method improves over the 
generalized gradient approximation^ (GGA) or lo- 
cal density approximation (LDA)— in the study of 
systems containing correlated electrons by introduc- 
ing the on-site Coulomb repulsion U. Results pre- 
sented here are consistent for values of U in the 
interval from about 3 to 5 eV, a reasonable range 
according to previous similar calculations of V^"*" 
compounds^"— and also of V'^^ compoundsi^ The 
calculations were fully converged with respect to 
the k-mesh and RmtKmoi- Values used for the k- 
mesh were 5x4x4 sampling of the full Brillouin zone. 
R.mtKmaa;= 6.0 was choscu for all the calculations. 
Selected muffin tin radii were the following: 1.72 a.u. 
for V, and 1.53 a.u. for O in the low-T structure and 
in the high-T phase, 1.77 a.u. for V, and 1.57 a.u. 
for O. 
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III. DISCUSSION OF THE RESULTS 

A. Structure 

We have taken the structural data for both the 
metallic phase (at 298 K) and the insulating one 
(at 120 K) from the work of Hodeau and Marezio^i^ 
The main structural features are depicted in Fig. 
[TJ The two phases comprise different V-V and V- 
O distances (the V-V distances are shown in Fig. 
[ijd), the V-0 distances are presented in Table 
As can be seen, in the V4-V2-V2-V4 chain of the 
metallic phase, the cations are equidistant. In the 
insulating state, the V2-V4 distance becomes shorter 
than the V2-V2 one, favoring a possible V2-V4 pair- 
ing. Even larger change occurs in the V3-V1-V1-V3 
chain distances from the high to the low tempera- 
ture structure, with the Vl-Vl distance becoming 
much shorter at lower temperatures. The distance 
between chains remains almost unchanged across the 
MIT. 

TABLE L Interatomic distances between V and O atoms 
(in A) in V4O7 for both the low-temperature (120 K) 
and high-temperature (298 K) phases. Structural data 
are taken from Ref. [T^ . 





120 K 


298 K 




120 K 


298 K 


d(Vl-02) 


1.81 


1.88 


d(V3-03) 


1.73 


1.78 


d(Vl-Ol) 


1.88 


1.89 


d(V3-06) 


1.92 


1.94 


d(Vl-Ol) 


1.92 


1.93 


d(V3-05) 


1.95 


1.96 


d(Vl-06) 


1.94 


2.05 


d(V3-07) 


1.99 


1.98 


d(Vl-04) 


2.05 


2.02 


d(V3-04) 


2.01 


2.03 


d(Vl-05) 


2.06 


2.03 


d(V3-05) 


2.15 


2.12 


d(V2-01) 


1.91 


1.90 


d(V4-06) 


1.90 


1.81 


d(V2-02) 


1.94 


1.95 


d(V4-05) 


1.96 


1.97 


d(V2-04) 


2.00 


2.03 


d(V4-07) 


1.98 


1.97 


d(V2-07) 


2.00 


2.02 


d(V4-02) 


2.00 


1.94 


d(V2-03) 


2.03 


1.97 


d(V4-07) 


2.08 


2.12 


d(V2-03) 


2.05 


2.01 


d(V4-04) 


2.13 


2.10 




FIG. 1: (Color online) (a) Pseudorutile blocks in V4O7. 
The lines represent the shear planes. We also show the 
magnetic structure of the low-temperature phase ob- 
tained for U= 4.1 eV (the coupling accross the shear 
planes becomes AFM for U > 5.4 eV). The ovals around 
VI pairs indicate spin-singlet states, (b) View of the 
chains where the zig-zag of the structure can be seen, 
(c) Unit cell of V4O7. (d) Table with the different V-V 
distances (in A) of the high-temperature (metallic) and 
low-temperature (insulating) structures at 298 K and 120 
K, respectively as obtained in Ref. I19I . 



the V'^^ and V'^ cations over the four cation sites 
and by the pairing of half of the V''"'" (those in VI 
positions) , and of all the V'^"'" in the short V-V bonds 
along the corresponding chains. The remaining un- 
paired V*"*" ions (V3), on the other hand, shift from 
the centers of the corresponding Og octahedra, form- 
ing one very short V-0 bond (V3-03) of less than 
1.8 A, as can be seen in Table HI 



According to the analysis performed in Ref. 19 by 
doing bond-length summations with the experimen- 
tal structural data, the high temperature metallic 
state is not comprised of completely disordered V'^"'', 
V^"*" cations. The 4224 chains contain around a 58% 
of the V'^"'" cations whereas the 3113 chains contain 
the same amount of the V^^ ones. Hence, in each 
chain, even in the metallic state, there remains some 
charge segregation, but a rather weak one. In the in- 
sulating phase, the 4224 chains become richer in 3-1- 
cations (about 83%), V2 having a valence of 3.24-I-, 
and V4 a valence of 3.09-I-. Meanwhile, the 3113 
chain has 83% of the V^^ cations, each of them hav- 
ing a valence of 3.83-I-. Thus, the low-temperature 
insulating state is characterized by the ordering of 



B. Electronic structure 

Bl. Low temperature phase 

The electronic structure of the low-temperature 
phase of V4O7 can be roughly understood by ana- 
lyzing the short-long V-V bonds alternating in the 
structure. These, in particular, modify the magnetic 
couplings, and also the distinct charge and orbital 
states for the V cations that such a structure might 
lead to. In order to analyze the electronic struc- 
ture of the compound, we have performed LDA+U 
calculations for several values of U (3.4 - 5.4 eV). 
For all the U values, the most stable solutions have 
in-chain antiferromagnetic (AFM) couplings. By in- 
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FIG. 2: (Color online) Magnetic couplings in the unit 
cell for the low-temperature structure. The AFM in- 
chain couplings and the inter-chain FM one (stable at 
lower values of U but very small in value) can be seen. 



creasing the U value, the inter-chain coupling, even 
though always almost negligible in value, switches 
from FM for lower U values to AFM for U > 5.4 
eV. The magnetic ground-state (consistent with the 
in-chain couplings for U ranging from 3.4 to 5.4 
eV) can be seen in Figs. [Tfa) and [2] Table HIl 
shows the magnetic moments of each V cation in the 
low- temperature structure obtained for U= 4.1 cV. 
From the charge/moment analysis, CO looks sub- 
stantial. In the V^+(3113) chains, the V magnetic 
moments inside the mufSn-tin spheres are about 0.8 
fiB, whereas for the V'^"'"(4224) chains the moments 
are of about 1.4 /is, slightly smaller than twice the 
value of V^^, that would be the case in a fully CO 
state. No noticeable difference in moment is ob- 
served between V2 and V4, in contrast to the result 
of bond length summations that predict a slight dif- 
ference in valence. The point charge model works 
very well, but hybridization also plays a role in de- 
termining the charge states and magnetic moments. 



TABLE II: Projection of the spin magnetic moments of V 
atoms inside the muffin-tin spheres in the V4O7 magnetic 
ground-state for the low temperature AFM insulating 
phase and for the high temperature FM metallic one. 



Low-T phase 




High-T phase 


Atom Magnetic Moment 


Atom 


Magnetic Moment 


Vla(b) 0.8(-0.8) fiB 


VI 


1.3 /is 


V2a(b) 1.4(-1.4) ^lB 


V2 


1.4 /is 


V3a(b) 0.8(-0.8) fiB 


V3 


1.2 flB 


V4a(b) 1.4(-1.4) hb 


V4 


1.2 /is 



Further indications of the large degree of CO, very 
close to a d^/d^ state in the corresponding chains, 
come from the study of the band-structure plots 
with band-character (the so called "fat-bands" plots) 
shown in the upper panels of Fig. |31 A band struc- 
ture plot is often a more reliable tool to analyze 
charge states in localized electron systems than a 



simple electron count inside the muffin-tin spheres4^ 
The total charge in both cases would be very simi- 
lar, causing the very small difference in total energy 
observed experimentally, but the electronic struc- 
ture and magnetic moments can be very different. 
This phase is an insulating one with a d-d gap of 
about 0.6 eV, for this particular value of U (4.1 eV). 
About 1 eV below the Fermi level, only d-bands are 
observed, with the O states much lower (starting 
at about -2.5 eV), as occurs in other vanadates like 
VO2,'''' SrVOar^ sympliflying the picture. We ob- 
serve, on the upper left two panels of Fig. [31 two 
cations with a d^ configuration (these are the cations 
in the V''~''-chain) and, on the right two panels, other 
two cations with a distinct d^ electronic configura- 
tion (those at the V'^+ chains). The "fat-bands" 
plots clearly show one band occupied for the V"*^ 
cations, and two for the V'^"'" ones. Close above the 
Fermi level, for each spin channel, six unoccupied 
bands can be observed, four of them from the V''"'": 
d^ cations and, split, two from the V'^"'": d^ cations. 

We can also describe the electronic structure of 
the material based on the partial density of states 
(DOS) plots of the various V cations in the low-T 
structure (see left panels of Fig. |4|). The differ- 
ent DOS plots for V^^ and V'^^ cations can also 
be observed, with a different position of the occu- 
pied bands with respect to the Fermi level: the rela- 
tively narrow V'*^ bands lie lower in energy, and the 
broader and more occupied bands corresponding to 
the cations are closer to the Fermi level. The 
DOS integrations inside the muffin-tin spheres give 
a less distinct picture of the CO phenomenon than 
the band structures, which are much clearer. 

In order to give further evidence of the large CO 
character of the low-temperature phase, we can con- 
sider a three-dimensional representation of the spin 
density; this also shows which orbitals are occupied 
at each site, i.e. it allows to find the 00. In Fig.jSfa) 
such representation is shown for U= 4.1 eV. 

Quite different spin densities (those of V'^+/V''+ 
cations) along the corresponding chains can be seen. 
The spin density of the V^+(d^) cations corresponds 
to that of a mixture of d^z and dyz orbitals, whereas 
for the V'*+(d^) it can be interpreted as a dxy orbital, 
using a local coordinate system with the z-axis di- 
rected from V towards its apical oxygen. 

In principle, one would expect the local octahe- 
dral environments of the V*"*" and V^"*" cations to be 
distorted accordingly, to accomodate one electron in 
the dxy singlet or two in the d^z, dyz doublet. How- 
ever, the distances we have summarized in Table |T] 
suggest complicated distortions of the VOe octahe- 
dra, which make it difficult to predict on the ba- 
sis of only local distortions which orbitals would be 
occupied )^^'^^ without analyzing the results of the 
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FIG. 3: (Color online) The upper panels show the band-structure plots with band-character of the majority states for 
the V cations in the low-temperature insulating phase, (a) Represents the bandstructure plot of the spin-up(down) 
channel of Vla(Vlb), (b) of the spin-up(down) channel of V3a(V3b), (c) of the spin-up(down) channel of V2a(V2b) 
and (d) of the spin-up(down) channel of V4a(V4b). The differences in the band-structure plots for the V'^^ (V2 and 
V4) and the V*+ (VI and V3) cations can be observed. The lower panels show the same plots for the majority (up) 
states for the V cations in the high-temperature metallic phase (a) for VI, (b) for V2, (c) for V3, and (d) for V4. 



calculations. 

On the basis of these results, one can suggest a 
possible explanation of the origin of the structural 
distortions found in V4O7 below Tc- ions oc- 

cupy two types of positions: VI in the middle of 
rutile-like chains, and V3 at the shear plane. When 
one electron occupies the dxy orbital in V^+ atoms, 
which occurs below Tc, there exists strong direct 
overlap of these orbitals along the chain, which can 
give rise to a metal-metal bonding of VI ions, like 
that in the singlet dimers in the insulating phase of 
VO24I However, the V^^ ions in V3 positions are in 
a different situation: they have a common face with 
V4 across the shear plane. In this situation, as ar- 
gued by Goodenough,^- the main instability of these 
ions would be the shift of such V^"*" from the center 
of its Og octahedron (ferro- or antiferroelectric dis- 
tortion). Indeed, as one sees from the Table HI just 
for this V there exists one very short V-0 distance, 
V3-03, of 1.73 A; all the other V-0 distances are 



larger and more or less equal. Thus, in the struc- 
ture of V4O7, with shear planes every four sites and 
with the obtained charge and 00, we would have an 
(anti)ferroelectric-type displacement of the "edge" 
V3 ions, but the "inner" V*^ ions VI would form 
short, presumably singlet (see below) dimers. The 
situation here resembles somewhat the M2 phase of 
V02,'^""^^ in which also half of V^+ ions form sin- 
glet dimers, but the other half has no dimerization, 
but has instead an antiferroelectric-type distortion 
("twisting") with the formation of one short V-0 
bond. 

The 00 in V'^^ ions, obtained in our calculations, 
is, however, quite different from that in V^^ cations: 
in V'^^ not d^j,, but rather two other orbitals, dxz 
and dyz, are occupied by two electrons in the low- 
temperature phase. These electrons do not have 
strong direct d-d overlap, but they can instead hop 
via bridging oxygens. By that process, the dxz elec- 
tron from one site hops to a dyz orbital of the neigh- 
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E(eV) E(eV) 

FIG. 4: (Color online) Left panels represent the par- 
tial DOS for each V cation in the low-T phase, where 
a charge disproportionation effect is clear. Right panels 
represent that of the V cations in the high-T structure; 
the different environment leads to a different DOS in 
each site, but the electron count is similar for all of them. 
Positive (negative) values indicate up-spin (down-spin) 
DOS. 



boring V, as sketched in Fig. [Bl where we show the 
V'^^ chain running along the c-direction in a more 
convenient projection. One sees indeed that e.g. the 
xz-electron from the site V4a can hop to a yz-state 
of the site V2b via a Pz orbital of the oxygen a, and 
then again to a;2:-orbital of the site V2a via the p^ or- 
bital of oxygen d. Similarly, there will exist hopping 

{yz)v4a - iPz)b - {xz)v2b - (Pz)c " {vz)v2a ■■■■ As a 

result, we would have two degenerate half-filled one- 
dimensional (ID) bands (containing two electrons 
per V^"*"). Such ID bands would naturally give rise 
to Peierls dimerization, which would open the gap in 
the spectrum and would decrease the energyJ^ This 
is a possible explanation of the formation of dimers 
(V4-V2) and (V2-V4) in V4O7. Similar dimeriza- 
tion was recently found in the rutile-like V^"*" chains 
in vanadium hoUandite K2V80i6,^'^ apparently due 
to the same mechanism. 



B2. Magnetic order 

To analyze the possibilility and the consequences 
of pairing of V atoms along each chain and to un- 
derstand the magnetic properties of the compound, 
we have calculated the strength of the magnetic cou- 
plings in the unit cell. Two kinds of couplings can be 




FIG. 5: (Color online) Three-dimensional plot of the spin 
density of the magnetic ground-states for both high-T 
and low-T phases of V4O7, isosurface at 0.1 e/A^ ob- 
tained using XCrysDen [5^. (a) Represents the plot for 
the low-temperature phase. Different colors (grayscales) 
represent spin-up(down) density. The upper panel in (a) 
shows the density plot for V and O atoms in the unit cell. 
The lower panel shows the CO character (d^/d^) along 
each chain in the structure (neither V nor O atoms have 
been depicted for simplicity), (b) Represents the spin- 
density plot for the high-temperature metallic phase. All 
V atoms exhibit a very similar spin density. 



considered in each chain (one short (s) and one long 
(1)) plus the inter-chain coupling (see Fig. [5]). We 
have used our total energy calculations for different 
coUinear magnetic configurations to quantify and de- 
scribe these couplings. To that end, we have fitted 
the total energies by a Heisenberg model, in the form 
H =\ Ji]SiSj, with S=l for the sites and 
S=l/2 for the ones. The sites would be 

expected to undergo singlet spin pairing whereas the 
V^"'"(d^) ones could pair in singlet, triplet or quintu- 
plet spin states. The values of the in-chain coupling 
constants, obtained in our calculations, are: J3; « 
40 K, Jss « 100 K, in « 10 K and J4S ~ 1000 K (3 
and 4 stand for V'^''" and V^^ chains, respectively, 
and "1/s" stand for the long and short bonds). 

The short-bond magnetic coupling between paired 
cations (J4s) is very large, w 1000 K, indicat- 
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FIG. 6: (Color online) Schematic representation of the 
possible hopping processes along the V"^"*" chains in the 
c-direction. Due to the electronic structure of these 
cations, the magnetic interactions are mediated by bridg- 
ing O p-orbitals, contrary to the direct V-V bond (not 
shown) that occurs in the V''+-V*+ chains. 



ing that most probably these ions form spin-singlet 
states, reminiscent of the ordering along the c-axis 
that happens in V02.-^' The situation is different 
in V'^"'" chains, where also the short-bond magnetic 
coupling {3 3s) is larger than that in the long bond 
(Ja;), but the difference is not as drastic as in the 
V'*"'' chains. This suggests that spin singlets would 
not be formed in the chains. The magnetic 
coupling strength between chains is quite small, but 
apparently sufficient to form long-range magnetic or- 
dering below Tat = 40 K. We can thus propose that 
the magnetic structure of this state may be the AFM 
ordering in V'^+ segments, coupled to ions V3 
with unpaired spins (see Figs. [2] and [Ha)). Due 
to the very small inter-chain magnetic coupling, two 
V''"'" cations per chain (V3) could be left partially 
disordered, even at low temperatures.^^ 

Previous works^L^ give further support to the im- 
portance of the pairing interactions as one of the 
driving forces for the MIT. In the work by Gossard 
et al-r^ the pairing in the chain (correspond- 
ing to our 34s) becomes complete only below 150 
K (100 K below the MIT transition). The magnetic 
coupling estimated using their nuclear-resonance 
frequency shifts is also large (J~ 500 K at 100 K). 
The nuclear magnetic resonance (NMR) spectrum 
of spin-echo amplitudes observed at 4.2 K implies a 
continuous singlet spin pairing in the AFM state, 
with the singlet V^^ spin-pairs unbroken by the 
AFM ordering. This is in agreement with our pic- 
ture. 



In addition, the possible formation of singlet spin 
pairs along the V'^+ chains is also discussed in Ref. 
|271 Due to the two echo lines appearing in the NMR 
spectrum one possibility is that i) one line arises 
from unpaired V'*^ with a moment of 1 hb, while 
the other comes from paired atoms which have 
reduced moments of 1 /i^j due to covalent singlet 
bonding of one of their two 3d electrons. The other 
option is: ii) both lines arise from unpaired A+ atoms 
that may have become crystallographically inequiv- 
alent in the AFM state. Structural results^^ at 120 
K, above the magnetic order temperature (Tjv), in- 
dicate the unpaired V^^ at the shear planes are not 
inequivalent. Heidemann et al^ gave futher con- 
firmation of the option i) and hence of the partial 
singlet spin pairing in the V^^ chains, possibly with 
some magnetism thereof remaining. 

From our calculations, we can argue that although 
Jss is larger than the coupling of V^"^ spins in the 
long bond, we do not expect real spin-singlets to be 
formed in the chains, but a more complicated 
spin structure is expected to develop. From our cal- 
culations, we can see that the magnetic moments 
obtained for the V^+(d^) atoms are reduced by hy- 
bridizations from the purely ionic picture. Also, the 
magnetic moment of these cations barely changes 
across the transition, as we saw in Table HIl Thus we 
propose that both V'^^ and half of V*"*" ions (those 
in V3 positions) participate in long-range magnetic 
ordering (see Fig. [TJa)). 



B3. High temperature phase 

The electronic structure of the high temper- 
ature/metallic phase does not present the large 
charge- ordering character, as the bond-length sum- 
mations correctly suggest We have performed 
LDA-t-U calculations for the same magnetic con- 
figurations used for the low-temperature structure 
with the same U values. The magnetic ground-state 
is now FM and all the sites are almost equivalent: 
the electronic structure is close to V^-^+ for all the 
atoms in the structure. The magnetic moments of 
the ground-state solution that can be seen in Table 
im are consistent with this picture. Some difference 
is seen in the magnetic moments due to the different 
band widths of the respective d-bands, caused by the 
different environments the V cations are in, but no 
trace of CO remains at high temperature, leading 
to a metallic solution. Also, by looking at the band 
structures with "fat-bands" (lower panels of Fig. [3]), 
we do not see the clear d^/(P ionic picture we ob- 
serve in the insulating low-temperature phase. The 
FM high-temperature band structure is consistent 
with the absence of charge-ordering we are describ- 
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ing, since it shows more evenly distributed occupa- 
tions in different bands. The spin density plot for 
the high temperature phase (see Fig. [SJb)) gives 
further evidence of this absence, confirming a non- 
integer charge occupation (roughly V'^'^+) at all the 
V sites. Even though a different DOS is seen for 
each of the V sites in the high temperature struc- 
ture (right panels of Fig. S]), the electron count is 
very similar for each of them, as the "fat bands" 
and magnetic moments also show. The FM order- 
ing (or paramagnetic state experimentally found at 
high temperatures) also increases the bandwidths 
when compared to the low-temperature CO AFM 
phase. The absence of CO is supported by experi- 
mental works. i^^i^ The high temperature suscepti- 
bility data for V4O , ^^1^^ are not conclusive about the 
sign of the Curie- Weiss temperature due to the high 
transition temperature. Our calculations suggest a 
change to a FM nature of the in-chain couplings in 
the high-temperature phase. One can qualitatively 
interpret the FM coupling in the high-temperature 
metallic phase as a double exchange interaction. 



C. Final remarks 

Thus we see that, in analogy with VO2, in the 
chains there occurs an 00 facilitating singlet pair 
formation. Also, the analysis of the electronic struc- 
ture helps to understand why only the inner (Vl-Vl) 
ions form such pairs, whereas (V4-V2) and (V2- 
V4) pairs are formed in V'^"'" chains. This has to do 
with the orbital occupation in the V'^+ chains, dif- 
ferent from that in the V^"*" chains: two orbitals, dxz 
and dyz, are occupied by two electrons at these sites, 
in contrast to dxy orbital occupation in V*"*" chains. 
Apparently, this is responsible for somewhat weaker 
dimerization (the changes in the V-V distance across 
the CO transition are, A{V4-V2) high-T,iow-T— 0.1 
A, as compared with A{Y1-Yl) high-T,iow-T= 0.12 
A, see Fig. Id). This also leads to smaller modula- 
tion of the exchange constants within the V'^"'" chains 
{1 short — 100 K vs 3 long = 40 K), whereas this dif- 
ference is much larger in V^"*" segments {Jshort = 
1000 K vs Jiong — 10 K). As a result, apparently 
V'^"'' ions do not form singlets, despite dimerization, 
but participate in magnetic ordering, together with 
unpaired V*"*" (V3) ions. This can explain the ap- 
pearance of magnetic ordering at a relatively high 
temperature Tn= 40 K. 



IV. SUMMARY 

In this paper, we have studied the electronic 
structure of V4O7 analyzing the differences be- 



tween the low-temperature (charge ordered, antifer- 
romagnetically ordered and insulating) and the high- 
temperature (charge disordered, paramagnetic and 
metallic) phases. 

Calculations show the strong CO in the low- 
temperature phase that, together with structural 
distortions, lead to a gap opening. The high- 
temperature structure is charge disordered, with all 
V atoms in a similar ~ d^'^ electronic configura- 
tion. All the in-chain magnetic couplings in the 
low-temperature phase are AFM, and the inter-chain 
couplings are very weak (and of different signs de- 
pending on the U- value) . Magnetic coupling changes 
to FM in the high-temperature phase as a result 
of the change in electron distribution caused by 
charge disorder. The AFM couplings in the low- 
temperature phase are obtained, being very strong 
in the V^+-V^+ short bonds, one order of magnitude 
stronger than the short V^+-V^+ magnetic couplings 
(and two orders of magnitude stronger than in the 
V'*+-V''+ long bonds). This suggests the formation 
of a spin-singlet between the "middle" Vs within the 
chains, leaving V3 atoms with unpaired spins. 
The couplings in the V'^"'' chains are more compli- 
cated, showing signs of slight moment reduction due 
to hybridizations. 

The 00 obtained in our electronic structure cal- 
culations (V+i dxy vs. V'^+: dxz Ay z) also helps to 
understand the different bond distance alternations 
that occur in the two different types of V chains: in 
the chains there is a strong metal-metal bond 
between the "middle" VI ions forming spin singlets, 
and an (anti)ferroelectric-type displacement of the 
"edge" V''"'" ions at the shear planes. On the other 
hand, the V^"*" chains, due to different orbital occu- 
pation, form two one-dimensional-like bands due to 
bridging with oxygens, that affects the type of cou- 
pling along these chains and leads, in particular, to 
a dimerization of V'^^ chains, which, however, is not 
strong enough to make the short V'^"'"-V'^"'" dimers 
spin singlets. 

As a result, most probably both the unpaired 
"edge" V^"*" and all V^"*" ions contribute to magnetic 
ordering below T^r , the detailed type of which is yet 
unknown. It would be very interesting to check the 
proposed magnetic ordering in V4O7 experimentally. 
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One should note that the one-electron picture of the 
electronic structure provided by DFT calculations 
cannot account for the formation of spin-singlets be- 
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